Pairing mechanism in multiband superconductors 
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It has been a long-standing puzzle why electrons form pairs in unconventional superconductors, 
where the mutual interactions are repulsive in nature. Here we find an analytic solution for renor- 
malization group analysis in multiband superconductors, which agrees with the numerical results 
exceedingly well. The analytic solution allows us to construct soluble effective theory and answers 
the pairing puzzle: electrons form pairs resonating between different bands to compensate the energy 
penalty for bring them together, just like the resonating chemical bonds in benzene. The analytic 
solutions allow us to explain the peculiar features of critical temperatures, spin fluctuations in un- 
conventional superconductors and can be generalized to cuprates where the notion of multiband is 
replaced by multipatch in momentum space. Therefore, finding effective attractions between elec- 
trons is no longer a necessity and the secret for higher superconducting temperatures lies in boosting 
pair hopping between different bands. 
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It has been 100 years since Heike Kamerlingh Onnes 
discovered the resistance of mercury suddenly drops to 
zero[T] when cooling down by liquid helium in 1911 and 
marked the birth of superconductivity. The exotic phe- 
nomena of superconductors remained mysterious until 
Bardeen, Cooper and SchriefFer (BCS)[J' came up with a 
complete theory in 1957. Despite the celebrating success 
of BCS theory, there are other superconductors remain 
queer [3] and cannot be explained solely by the electron- 
phonon interactions, including cuprates [H [S], heavey- 
fermion compounds [MS] , organic superconductors and 
recently found iron-based superconductors fO] [1^. Per- 
haps the most intriguing puzzle for these unconventional 
superconductors is the pairing mechanism: what is the 
glue to pair up electrons from mutual repulsive inter- 
actions? The experimental evidences suggest that pair- 
ing in the unconventional superconductors is not due 
to electron-phonon interactions. Due to strong mag- 
netic correlations [TTHIl] in these materials, it is pro- 
posed that spin fluctuations 1 1 SffTO] may play the role of 
glue to pair electrons up. Recent renormalization-group 
(RG) studies[22] indeed reveals the close relation between 
spin fluctuations and unconventional superconductivity 
in iron-based materials. 

In this Letter, we investigate the pairing mechanism 
in multiband superconductors by the unbiased RG ap- 
proach. In general, the RG equations are coupled non- 
linear first-order differential equations and make any sim- 
ple understanding beneath the messy numerics inacces- 
sible. However, making use of classification scheme by 
RG exponents, we show that the dominant interactions 
are intra-band g and inter-band g^_ Cooper scattering. 
It is rather surprising that these two dominant couplings 
are captured by a set of analytic solutions. The solu- 
tions are elegant and simple enough to reveal the pairing 
mechanism in multi band superconductors. 

The binding energy for Cooper pair formation is dic- 
tated by a small parameter 5 = g — \g±]. As long as the 



inter-band pair hopping is larger than the intra-band. 
Cooper pairs resonating between different bands become 
stable despite of the repulsive intra-band interaction. 
The picture of resonating Cooper pairs between differ- 
ent bands leads to several natural consequences. First of 
all, though the coupling strengths of g and g^_ are large, 
the binding energy of Cooper pairs (and thus the critical 
temperature) is determined by their difference 5, which 
is one order smaller than the bare couplings as detailed 
later. Secondly, it can be shown that the inter-band pair 
hopping also give rise to spin fluctuations at the nesting 
vector which connects the dominant Fermi surfaces. The 
inter-band pair hopping not only explains why strong 
magnetic fluctuations often show up in unconventional 
superconductors but also pins down at what momentum 
the spin fluctuations should appear. 

To see how the analytic solutions emerge from the 
RG analysis, we choose the iron-based superconductor as 
an demonstrating example. We start with a five-orbital 
tight-binding model pPH^ for iron-based superconduc- 
tors with generalized on-site interactions, 
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where 



1,2, 



, 5 label the five d-orbitals of Fe, 
1 : d^z^-B.^, 2 : dxz, 3 : dyz, 4 : dx^-Y^i 5 : dxY- 
The kinetic matrix Kab in the momentum space has been 
constructed in previous studies[3D]. The generalized on- 
site interactions consist of three parts: intra-orbital Ui, 
inter-orbital U2 and Hund's coupling Jh- To simplify 
the RG analysis, we sample each pocket with one pair of 
Fermi points (required by time-reversal symmetry). This 
is equivalent to a four-leg ladder geometry with quan- 
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FIG. 1: (a) Fermi surfaces for the five-orbital Hamiltonian at 
electron doping x = Q.\ with ladder geometry. The RG expo- 
nents for (b) Cooper scattering (c) forward scattering in spin 
sector (d) forward scattering in charge sector are presented. 
The dominant bands are M,N — 1,2, while the subdomi- 
nant ones are m,n = 3,4,5. The interaction profile is set to 
Ui = 2.8 eV, U2 = 1.4 eV and Jh = 0.7 eV. 



tized momenta and the effective Hamiltonian consists of 
five pairs of chiral fermions as shown in Fig. 1. 

In weak couDling |23H25j . the allowed interactions are 



Cooper scattering cf^ , cf ■ and forward scattering /^^ , f[j 
between different bands, where u, p denote the spin and 
charge channels respectively. We integrate the coupled 
RG equations numerically and find all couplings are cap- 
tured by the scaling Ansatz p51l27j . 
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where Gi is an order one constant and < 7i < 1. The 
scaling exponent 7^ help us to build the hierarchy of all 
relevant couplings without ambiguity. 

In the doping range x = to a; = 0.12, these exponents 
are shown in Fig. 1. The couplings with 7^ = 1 are the 
most dominant, including intra-band Cooper scattering 
c^i , C22 < and the inter-band ones c^2 > between 
the hole pocket centered at (0, 0) (band 1) and the elec- 
tron pocket at (±7r, 0) (band 2). The positive sign of the 
inter-band Cooper scattering C12 leads to the s± pair- 
ing symmetry. Note that our RG Ansatz predicts the 
dominant superconducting bands with the correct pair- 
ing symmetry as obtained by the fRG method. In fact, 
a detail analysis including all subdominant relevant cou- 
plings lead to correct signs for all gap functions in differ- 
ent bands. 

Introduce the couplings, c = {cii+ C22+ c1i + (^22) I ^ fo'" 
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FIG. 2: Comparison between the analytic solutions for g{l) 
and g^{l) and the numerical RG flows for c{l) and c^{l). 



intraband pair hopping and c± = (c^2 + c: 



-12' 



-c?i)/4 



for interband pair hopping. It is remarkable that the RG 
flows obtained in numerics, as shown in Fig. 2, are well 
captured by the analytic solutions g{l) and g±{l), 
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where l± = — 1/ [17(0) ± .gj^(O)] can be extracted nu- 
merically. Because the RG flows diverge a.t I = 1,1, it 
flxes I- = Id- Meanwhile, it is reasonable to require 
that c(0) -|- c_l(0) = g{0) + g±{0), ensuring the coupling 
strengths are of the same order. Thus, the two parame- 
ters l± can be determined without any fitting. 

What is the secret message behind the analytic so- 
lutions? Consider a two-band BCS Hamiltonian with 
intra-band g and inter-band g± pair hopping. Though 
not widely known, RG equations for these couplings can 
be derived from the dependence of the gap functions. Af- 
ter some algebra, the RG equations for g and g± are 



(4) 



The above non-linear coupled equations can be solved 
exactly, giving the analytic solutions we discussed pre- 
viously. Therefore, the effective theory for iron-based 
superconductors is the multiband BCS Hamiltonian, 
proven by matching the RG flows together. 

RG flows for g and g± are shown in Fig. 3. If the 
intra-band pair hopping is larger than the inter-band 
one {g — g± > 0), the couplings flow toward the Fermi- 
liquid fixed point. If the inter-band pair hopping is larger 
{g~9i- < 0)1 they flow toward the superconducting phase 
with s± pairing symmetry. It is worth mentioning that, 
if the initial couplings are close to the symmetric ray 
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FIG. 3: (a) RG flows for two-band BCS Hamiltonian hosting 
three difi'erent phases: Fermi liquid, conventional and uncon- 
ventional s±-wave superconductors, (b) Typical RG flows for 
bare couplings close to the symmetric ray g± = g. The cor- 
responding trajectory on the g — gi_ plane is also illustrated 
in part (a). 



g ~ g±, it will flow toward the Fermi-liquid fixed point 
first and then turns around to the unconventional super- 
conducting state. This means that, upon cooling down 
the system toward the critical temperature, it exhibits 
non-trivial crossover properties over a wide range of tem- 
peratures as shown in Fig. 3(b). However, in the absence 
of inter-band pair hopping {g± = 0), a negative ^(0), i.e. 
attractive interaction, is required to trigger the super- 
conducting instability. The RG flows in the special case 
g± = give rise to the commonly accepted criterion that 
an effective attraction is necessary for Cooper pair for- 
mation. However, the criterion is obviously wrong as the 
RG flows in the upper plane is quite different from those 
in the horizontal axis. 

The key parameter for Cooper pair formation in multi- 
band superconductors is 6 — g — \g±\. Despite of the 
energy penalty g to form a Cooper pair within the same 
band, through inter-band pair hopping, a Cooper pair 
gains — |5_l| benefit through resonating between different 
bands, just like the resonating chemical bonds in ben- 
zene. Thus, attractive interactions are no longer nec- 
essary. Following the textbook calculations, the critical 
temperature is 



/crT^- 1.14 Ae-^/l'^l, 



(5) 
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where A is of the same order of electronic band width. 
Setting A = t = leV, reasonable estimate for the hopping 



FIG. 4: Critical temperature Tc versus the interaction 
strength. The on-site interaction strength is Ut = U1+U2 + JH 
with the interaction proflle U1/U2 = 2 and Ui/Jh ~ 4. The 
critical temperature is dictated by the small parameter 5, 
which can be extracted from numerical RG flows and is one 
order smaller than the bare interaction strength Ut- 



amplitude i, we extract d from the numerical RG flows 
for different on-site interaction strengths Ut = Ui + U2 + 
Jh and plot the critical temperatures in Fig. 4. The 
parameter 5 is one order of magnitude smaller that the 
bare coupling Ut/t. If one chooses Ui = 2.8, C/2 = 1-4 
and Jh = 0.7, the predicted critical temperature is about 
56 K - quite a reasonable estimate. 

The inter-band pair hopping g±^ brings in another in- 
teresting property in unconventional superconductors. 
Making use of the instability analysis developed by Wang 
and Lee 22 , the interband pair hopping also enhances 
spin-density-wave (SDW) instability if the momentum 
Q — K1 — K2, connecting the two Fermi surfaces, is close 
to half of the reciprocal lattice vectors, i.e. Q — G/2. In 
the case studied here, Q — (tt, 0) satisfies the above con- 
dition. Therefore, it is expected that the antiferromag- 
netic spin fluctuations at the momentum Q are enhanced 
along with unconventional superconductivity. 

Generalizing Shankar's seminal work^U] in two di- 
mensions, we consider a two-pocket model with generic 
4-fermion interactions including intra-band forward 
scattering -Fpp(0i,02), inter-band forward scattering 
Fpp{9i,62), intra-band Cooper scattering Cpp{9i,02) 
and inter-band Cooper scattering Cpp{9i,92), where 
P ~ 1,2 is the band index with the convention (1,2) = 
(2, 1) and 9i represents the angle of the corresponding 
momentum. Detail derivations of the RG equations will 
be given elsewhere. In the absence of inter-band pair 
hopping, we have checked that our RG equations reduce 
to those derived by Shankar. Under RG transformation, 
the forward scattering does not flow but the Cooper scat- 
tering is described by the RG equations. 



Cpp{9i,62) — 
Cpp (6*1,02) = 



d9 

dd 

2^ 



[Cpp{0i,e)Cpp{e, 9-2) + Cpp{ei,9)Cpp{e, 62)] , 
[Cpp{ei,9)Cpp{9, 02) + c'pp{ei,0)Cpp{e, 02)] 



(6) 



Assuming the density of states and the bare couphngs 
are rotationally invariant, i.e. G{9i,02) = G{0i — 02), 
the RG equations can be decoupled by the partial-wave 
decomposition. By identifying Cpp — >■ g and Cpp — > g± , 
the same set of RG equations as in Eq. HI appears and 
leads to similar unconventional superconductivity in two 
dimensions. 

If both Fermi surfaces in the two-pocket model locate 
at Kj = 0, the superconductivity is still driven by the 
inter-band pair hopping but it does not come with any 
enhanced antiferromagnetic spin fluctuations. To some 
extent, one can say the spin fluctuations are "symptoms" 
of unconventional superconductivity but not cause of it. 
It would be truly exciting to look for realistic multiband 
superconductors with the band structure discussed here: 
unconventional superconductivity without enhanced an- 
tiferromagnetic spin fluctuations. 

Our previous calculations concentrate on the dominant 
two bands. But, it is reasonable to include all active 
bands in the generalized BCS Hamiltonian when quanti- 
tative accuracy is required. In addition, for single-band 
materials but with significant variations in density of 
states at different momenta, the pairing mechanism dis- 
cussed here is also at work. We have carried out primi- 
tive RG analysis on cuprates by cutting the single-band 
Fermi surface into 16 patches (minimum to differenti- 
ate nodal and antinodal regimes with four-fold dihedral 
symmetry). The dominant patches locate at (7r,0) and 
(0, tt) with sign-reversed gap functions, agreeing with the 
d-wave symmetry and the enhanced spin fluctuations at 
(tTjTt). How to rigorously generalize the analysis pre- 
sented here from multiband superconductors to multi- 
patch single-band system remains open. Though the an- 
alytic solutions show that the couplings flow toward the 
multiband BCS Hamiltonian, which can be treated easily 
at mean-field level, the strong-coupling theory may not 
be of the same origin. However, the elegant and simple 
analytic solutions extracted here may provide some help- 
ful hints in constructing the appropriate model or even 
the ground-state wave function. 
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